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The zero-temperature Casimir-Lifshitz force between two plates moving parallel to each other at arbitrary 
constant speed was found in [New J. Phys. 11, 033035 (2009)]. The solution is here generalized to the case 
where the plates are at different temperatures. The Casimir-Lifshitz force is obtained by calculating the elec- 
tromagnetic stress tensor, using the method employed by Antezza et al. [Phys. Rev. A 77, 022901 (2008)] for 
non-moving plates at different temperatures. The perpendicular force on the plates has contributions from the 
quantum vacuum and from the thermal radiation; both of these contributions are influenced by the motion. In 
addition to the perpendicular force, thermal radiation from the moving plates gives rise to a lateral component 
of the Casimir-Lifshitz force, an effect with no quantum-vacuum contribution. The zero-temperature results are 
reproduced, in particular the non-existence of a quantum-vacuum friction between the plates. 

PACS numbers: 42.50.Nn, 42.50.Lc, 46.55.+d 



I. INTRODUCTION 

The theory of the electromagnetic force between macroscopic bodies due to vacuum zero-point energy was developed by 
Lifshitz and co-workers ||Tl|2l[3i, following the consideration by Casimir of perfect mirrors [4 1 . Recently, Lifshitz's formula IHO 
for the quantum-vacuum force between two parallel dielectric plates was generalized to the case where the plates have an 
arbitrary constant lateral motion [5|. The motion introduces considerable extra complexity into the calculation but the final 
answer turned out to be relatively simple (it will be recovered here as a special case). As described in |5|, a contradictory 
literature had built up around this problem, with claims of a lateral ("quantum friction") force on the plates but disagreement as 
to its magnitude. The exact solution |5 1 showed that, whereas the perpendicular force is altered by the motion, there is in fact no 
lateral force on the plates. 

Lifshitz ||T][3l also included the effect of thermal radiation in his analysis. The Casimir-Lifshitz effect is therefore also taken to 
describe forces that have a contribution from thermal radiation as well as from the quantum vacuum. The formalism developed 
by Lifshitz, however, cannot be used for plates at different temperatures. The problem of non-moving parallel plates at different 
temperatures has been comprehensively studied by Antezza et al. ||6|. These authors used Rytov's theory of electromagnetic 
fluctuations 1 7 1 to compute separately the contribution of each plate to the electromagnetic stress tensor, which gives the Casimir- 
Lifshitz force. The approach is similar to Polder and van Hove's |8| treatment of radiative heat transfer between plates. The 
Casimir-Lifshitz effect out of thermal equilibrium has been investigated experimentally by measuring the surface-atom force 
between a Bose-Einstein condensate and a dielectric substrate [9J; the measurement was done for a substrate temperature higher 
than that of the surrounding environment and results were in agreement with theory ifTOll . 

In this paper we use the method of Antezza et al. {E\ to compute the Casimir-Lifshitz force between plates at different 
temperatures when one of the plates moves at an arbitrary constant speed parallel to the other The zero-temperature result [Jl 
therefore emerges as a special case. As a consequence of the thermal radiation, the electromagnetic stress tensor has non-zero 
off-diagonal components which cause a lateral force on the plates. The lateral force acts against the relative motion and is 
present even if the plates have the same (non-zero) temperature. Thus, the thermal radiation provides the means of restoring 
thermodynamic equilibrium, which requires the damping of the relative motion. The thermal radiation from each plate has a 
Planck spectrum \/[exY>{huj /kgT) — 1] in a frame co-moving with that plate and since the frequency changes under a Lorentz 
transformation, the spectrum of the radiation from each plate is different in frames with relative motion. This change in the 
thermal spectrum of the radiation from one plate when viewed in a frame co-moving with the other causes the lateral force 
between plates at the same temperature. The disappearance of the lateral force at zero temperature can therefore be attributed to 
the Lorentz invariance of quantum- vacuum energy 

The purpose of this paper is to derive the exact Casimir-Lifshitz force for this problem; the plates are at different temperatures 
and have arbitrary electric permittivities and magnetic permeabilities, and arbitrary constant lateral motion. We will not here 
enter into considerations of specific materials or approximations that hold in particular regimes, since the exact solution already 
requires a lengthy derivation. Moreover, the subject of Casimir-Lifshitz forces on moving media has suffered from a conflict of 
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FIG. 1 : Two plates at different temperatures, with general electric permitivities and magnetic permeabilities, lie in the yz-plane with constant 
separation a. The second plate moves in the postive y-direction with speed /3c, where c is the speed of light. We compute the Casimir-Lifshitz 
force on Plate 1 . 



assertions based on various attempts at approximations (see references and discussion in fSl) and we wish to make as clear as 
possible that the problem considered here is treated in full generality, within the framework of macroscopic electromagnetism. 
The geometry and notation of the problem are set out in Figure [T] and caption. The plates have are taken to have a separation 
a in the x-direction. Plate 2 moves in the positive y-direction and we use /3 to denote its speed in units of the speed of light c, 
i.e. the speed is (3c. The formalism used in the calculation is described in Section [IH it is essentially that of Antezza et al. |6|, 
except that we allow for a magnetic response in the plates. This results in an expression for the contribution of each plate to 
the electromagnetic stress tensor in terms of a Green tensor of the vector potential. The contribution of the non-moving plate is 
found in Section |inj this involves the computation of the required Green tensor In Section |lVj the contribution of the moving 
plate is found from that of the non-moving one by taking account of its differing material properties and its motion. Finally, the 
Casimir-Lifshitz force is obtained in SectionlV] 



II. FORMALISM 



Casimir-Lifshitz forces are given by the expectation value of the electromagnetic stress tensor |[Tl|2l[3 

(T = (D ® E) + (H B) - il((D • E) + (H • B)). (1) 

For zero-temperature problems the averaging in ([TJ is simply the (electromagnetic) vacuum expectation value, whereas if thermal 
effects are included it also contains the thermal averaging. Since the plates are at different temperatures the contribution of each 
to the fluctuating electromagnetic fields is computed separately. The relevant components of the stress tensor are a^x^ which 
gives the perpendicular force on the plates, and a^y {— Cyx) which gives any non-vanishing lateral force on the plates. We 
denote the contributions of Plate 1 at temperature Ti to these components by <Jixx{Ti) and (yixy{Ti), with a similar notation for 
the contributions of Plate 2 at temperature T2. Thus, we write the relevant components of the stress tensor as 

(^xx{Ti,T2) — uixx{Ti) + a2xx{T2), (2) 

T2) = <Jlxy{Ti) + (T2xy{T2). (3) 

We use a gauge in which the scalar potential is zero, so that the electric and magnetic fields are given by 

E = -(9f A, B = V X A, (4) 

and the vector-potential wave equation is 



A(r,w) = / dtA(r,t)e"^', (5) 

J — OO 

— -V X -eH^ 



1 w^' 

V X V X ~e{uj) ^ ) A = ^oj, (6) 



3 



where j is the current density of any sources. The retarded Green tensor for the vector potential satisfies 



1 ur 
V X — -V X -e(uj)^ ) G(r,r',tj) = l(5(r-r'). 



and relates the vector potential at r to the current density at r': 

A(r,c^)=/io / dVG(r,r',c^).j(r',c^). 



(7) 



(8) 

In Rytov's formaUsm Q, the fluctuating electromagnetic fields are regarded as being produced by fictitious fluctuating currents 
through ([8]). 

The expectation values that occur in the stress tensor ([T]) can all be calculated from an equal-time correlation function of the 
vector potential operator To avoid ambiguities in the operator ordering we follow the usual quantum prescription and consider 
the correlation function 



Using (jSj) and ([8]), we can relate (j9]l to a similar correlation function of the fictitious currents: 



(9) 



(ia(r,t)i6(r',i))s =■ 



Mo 



dJ / dV / d3r"'Ga,(r,r",a;)Gbrf(r',r"',c.') 



(2^)2 

X (j,(r",c.)jrf(r"',c.'))5e-'('^+-')*. 
The correlation function of the currents is given by the fluctuation-dissipation theorem 13] fTT) as 



(10) 



(j,(r",c.)3,(r"',a;'))s = 



Mo 

X (5(w + J) coth 



Im(^(cj)) 



2 Im(£(o;))(5,, + ^-^i^P^ (V" • V"'<5,rf - V','V7) 



(5(r" - r'") 



2kBT 



(11) 



The right-hand side of ( 11 1 is constructed from the differential operator in (j7|l that functions as the "inverse" of the Green tensor, 
which in turn serves as a generalized susceptibility ll3]|TT] for the vector potential through ([8]l. The derivatives in ( 11 1 act on the 
delta function S{r" — r'") and have been written in a form that makes clearer the symmetry under interchange of r" and r'"; 
throughout this paper we use Cartesian coordinates, for which V"(5(r" — r'") = — V"'(5(r" — r'"). Equation (11 1 is usually 
written without the term containing the magnetic permeability, as in |61 for example. Since we allow for arbitrary magnetic 
permeabilities in the plates we must write the more cumbersome expression ( fTT) . The final formulae for the contributions to 
the electromagnetic stress tensor will all be written in terms of reflection coefficients for transverse electric (TE) and transverse 
magnetic (TM) plane waves, and the permittivities and permeabilities will appear only through these reflection coefficients. 
Not surprisingly, the same result is obtained if one ignores the magnetic response in the derivation but then includes it in the 
reflections coefficients once the final expression is obtained. 
Using ( [TT] i and the Green tensor property 1 3 , 5 1 



Gm{v\ r'", -Lo) = GUv', r'", c^) (to real) 



{Aa{r,t)Ai,{r',t))s= / dtu {AA)ab{r,r' ,uj), 



we can write ([TO]l as 



where {AA)ab{j^, r', w) is a frequency-space correlation function defined by 



(12) 



(13) 



(AA)af,(r,r',w) 



27r 



d^r" / d^r'" coth 



\2kBT 



-^Im(£(cj))(5cd + I , (V -V dcd~V^V^ 



Gac{vy,Uj)GUv'y\u:) 

5{v" -v'"). 



(14) 
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(We note in passing that the familiar property of the Green tensor for non-moving media, namely that G6d(r', r'", w) is equal to 
Gdb{^"' ; r', does not hold for the problem considered here since it requires time-reversal invariance |3||.) E quation ( 12 1 and 
the similar property of the permittivity and permeability ||T4| allows us replace duj by 2Re dw in (|13|): 

POO 

(ia(r,t)i6(r',i))s = 2Re / c^lj (AA),b(r, r', w), 



(15) 

Correlation functions for the electric and magnetic fields are easily obtained by using (j4]) and (jSjl in the foregoing derivation, 
with the result: 

/■oo 

(16) 
(17) 

where Cacd is (in Cartesian coordinates) the permutation symbol. Equations ([T6]l-([T7]i are suitable for finding the stress tensor 
([T]l between the plates, where D — SqEj and H = B//io. 

The calculation proceeds by using (14 1 to calculate the contribution of each plate to the correlations functions (16i-(17i 
between the plates, and thence to the stress tensor components (j2]) and (j3]l. Everything is determined once we find the relevant 
Green tensor in ( 14 1. We first consider the contribution of Plate 1, for which the source points r", r'" in ( 14 1 are located inside 
Plate 1 while the points r, r' are placed between the plates. 



Jo 

/•oo 

{Ba{r,t)Bb{v',t))s = 2Re / duj eacdebef^cK{^^)df{r,r' ,lu), 



m. CONTRIBUTION OF PLATE 1 



A. Green tensor 



The required Green tensor can be written down using a method previously employed to solve the zero-temperature case ||5|, 
and also to find the Casimir-Lifshitz forces between birefringent plates (TT\. The solution of the monochromatic equation (j7| has 
a simple physical meaning: an oscillating dipole at the point r' emits electromagnetic waves of frequency uj and G(r, r', uj) is 
the resulting vector potential at the point r. The second index in Gab represents the orientation of the dipole at r', while the first 



index represents the components of the vector potential at r. Here we require the Green tensor G(r, r" ,lu) in ( 14 1, where the 
dipole sits at r" inside Plate 1 and the resulting vector potential is measured at r between the plates. To write down the solution 
we exploit the homogeneity of the problem in the y- and z- directions to Fourier transform the Green tensor as follows: 

/•oo /'OO 

G{x,x",u,v,io)= dy dzG(r,r",w)e-™(f-f")-™("-"") (18) 



so that we decompose the waves emitted by the dipole into plane waves. In the absence of the plate boundary, the solution of (j?]) 
inside the material of Plate 1 is the bare Green tensor |3 1 (with source point labeled r") 

, ' (19) 

^-tw{x-x )g_^ 2; < x" 

^ f k± ® k± + £1 (c^) ^ 1 V (20) 



= (±wi,w,i;), = Wei(w)/ii(w)^ - s2, s = \J u?- + v'^ . (21) 



The two possibilities in ( 19 1 represent plane waves propagating to the right (first line) or to the left (second line), with wave 
vectors k* . When x" is in Plate 1 and x is in the gap between the plates, only the right-moving plane wave can reach x from 
x" and it does this by first passing though the plate boundary into the gap. Once in the gap the plane wave can propagate to x 
directly, or after an even or odd number of reflections off the plates, giving right- and left-moving plane waves in the gap with 
wave vectors 



/ up- 

(±w,w, w), -s2. (22) 
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The solution for the Green tensor is the Hnear superposition of all these possibiUties; by the same reasoning employed in ||5] and 
fl2J we can write it as 



Gf II \ { iwx I 

(x,x = ( e + 



-iw{x — ^i) iwa-ry 



1 



(23) 



In ( |23] l, T is the transition operator (matrix) for the right-moving plane wave in Plate 1 as it passes into the gap, R2 is the 
reflection operator for a right-moving plane wave at plate 2 and Ri is the reflection operator for a left-moving plane wave at 
plate 1 . The inverse matrix in ( 23 1 is a geometric series representing every possible number (including zero) of double reflections 
off both plates, the exponentials providing the propagation distance 2a for each double reflection. Each term in (23 1, after the 
series expansion of the inverse matrix, has an overall exponential factor that accounts for the propagation distance involved. 

The operators T, R2 and Ri were derived in |5| by decomposing the plane waves into TE and TM polarizations in a frame 
co-moving with the relevant plate; they can be written 



tsriBi (g) riBi 



Ri = TEinEi (gi nEi + rsinj^i n^^, R2 = rE2nE2 ^e2 + '^32^^^ ® ^32- 



(24) 
(25) 



The scalar coefficients in ( 24 1 and ( 25 1 are transmission and reflection coefficients for TE (E) and TM (B) polarized plane 
waves 1 13 1 evaluated in a frame co-moving with the relevant plate ||5l: 



tE 



2wi 



= 1 - tei, 



tB 



TEl 



rE2 



fi2{u;')w - W2 
IJ.2{i^')w + W2 



rB2 



fj,i{uj)'w ~ wi Si{lu)w — Wi 

fJ,l{Lu)'W + Wl' ei{Lu)w + 'Wi 

e2iuj')w ~ W2 



(26) 
(27) 

(28) 
(29) 

Note that the reflection coefficients ( |28] l at Plate 2 are those in a frame co-moving with that plate; hence the Lorentz-transformed 
frequency ui' and y-component of the wave vector u'. The unit vectors in in (24 1 and (25 1 are given by 13] 



W2 = \ £2{uj')^i2i^')- 



e2{oj')w + W2 ' 

u;' = 7(a;-/?cu), u' = 7(1* - ^w/c), 7 = (1 - /J^)- 




n_Bml 




v{—uj/c + Pu) 



n 



B2 



±{~s'^ + Puju/c) 
w(u — f3uj/c) 
wv 



(30) 



(31) 



(32) 



CJ^.s2-2/3wu/c-/32(t;2-cj2/c2) 

and are polarization directions for plane waves, where the superscript + (— ) refers to a plane wave moving to the right (left). 
The Green tensor (23 1 is now completely specified, but before using it to compute the correlation function ( 14i we simplify the 
inverse matrix in ( |23| l. This is done in the Appendix, using the procedure employed in |5l- 



B. Contribution to the stress tensor 



In ( fTS) and ( [23] l one can see the entire spatial dependence of the Green tensor; spatial derivatives of the Green tensor result in 
factors of wave-vector components. We insert the fourier-transformed Green tensor ( 18 1 in ( 14 1 and integrate by parts so that the 
derivatives in (14 1 act on the Green tensors; these derivatives of the Green tensor are evaluated by multiplying by the relevant 
wave-vector component. The integrations with respect to r" and r'" are over the volume of Plate 1 ; we drop the boundary term 
at the surface of the plate that results from the integrations by parts, since we can take the fictitious currents to vanish on the 
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surface without affecting the contribution of the plate. This resuhs in the following formula for the correlation function ( [T4) l: 



(AA)ab(r,r',cj) 

47^3 



dx" 



du I dv coth 

-oo 



Jli(y-y')+iv(z-z') 



Im(^i(a;)) 



X Gac{x,x",u,v,uj)Gla{x',x",u,v,uj) <! ^lm(ei(w))(5cd + + s'^)Scd 

-|wip(5rf2;(5c2; - U^SdyScy " W^'^dz'^cz " Wiuddx^cy - wluSdyScx - WivSdxScz - wlvSdzSc 
-Uv{SdyScz + SdzScy)] 



(33) 



The expression for the Green tensor given in the Appendix is now used in ( 33 1 to compute the correlation function. As is clear 
from the simple dependence of the Green tensor on x", the integration with respect to x" in ( 33 i is performed using 







r dx" 




' — oo 





1 



21m(wi)' 



(34) 



Electric and magnetic field correlation functions follow from (16 1 and (17 1, and with r' ^ r these determine the contribution to 
the stress tensor (jT|) (where D = EqE and H — B//io between the plates). The fact that the real part is taken in ( 16 1 and ( 17 1 
means that the result can be split into terms from evanescent waves (a; < cs) with imaginary w and propagating waves (w > cs) 
with real w (see (|22]l). Use of the identities 



Im(ei) = 2Re{^i)Re{wi)lm{wi) — lm(/ii) 



[(ReK))2 - (Im(u;i))2] + 



Re{el^ilwi) = -^Re(wi) + s^) , 



(35) 
(36) 



c 

lm{£lnlwi) = Im(u;i) (| 



(37) 



allows the explicit occurrence of wi, e\ and to be re-expressed in terms of the reflection coefficients. With the definitions 



aEE ■■= 1 - e^'-'-rBirss, ass := 1 - e'"''^ r bxt bi. 
aEB 1 - e2»-r£irB2, asE 1 - e^^'^Vsir-^s, 



the result for (Jixx{Ti) and (Tixy{Ti) can be written 



O'lxxiTi) 

h 



167r3 



du I dv j du coth — — — 
-oo Jo \2kBli 



X { 4Im(r£;i) — ^ '-^ ^ 2 h (b ^ s) 

y \{cs'^ - u(3uj)'^aEEaBB + aEBaBEv'^fi'^c^w'^l 



167r3 



du I dv I dio coth , 

-oo J cs \2kBll 



w [{cs^ - uf3ujf + v'^0^c^w^ 



X { [^-\rEi\ 



(l + \rE2\^) \aBB\Hcs' - uPuf + (l + |rB2|') Wbe^v^P^ 

|(cs2 - u(3ujYaEEaBB + aEBO-BEV^P^C^W^ 



(38) 
(39) 



(e^b)}, 



(40) 
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O'lxyiTl) 



t poo poo pes / t 

' du I dv I diu coth ( J^]e-^-\^\u lies''- uPio)^ + v^P^c'w^ 



^ poo I^OQ /'OO 



X U 1 - kisi 



(l - |ri52|') |asBp(cs2 - u/3u;)2 + (l - |rB2|') \aBE\'v'P' 
|(cs2 - uPuYaEEa-BB + aEBaBEv'P'c'^w^l^ 



+ (e ^ b) 



(41) 



The first integrals (|40]) and ( |4T] i come from the evanescent waves while the second integrals are from the propagating waves. It 
is very important to note that, in addition to their explicit appearance in (40 1 and (41 1, the reflection coefficients at each plate are 



present in the quantities ( 38 1-( 39 1 



It is clear on symmetry grounds that the contribution (40 1 to the perpendicular force on the plates cannot depend on the sign 



of p. One easily sees this to be the case in (40 1 by considering a Taylor expansion of the integrands in powers of Odd powers 
of (3 are accompanied by odd powers of u and so these terms vanish after the integration with respect to u; only the even powers 
of /3 in the expansion, which are accompanied by even powers of u, contribute to aixx{Ti). (Note that the reflection coefficients 



(28 i-(29 1 at Plate 2 must also be expanded in powers of (3 to obtain the Taylor expansion of the integrands.) Similar symmetry 



considerations show that the contribution (41 1 to the lateral force must change sign when (3 changes sign. This is also seen to be 



the case by a Taylor expansion of the integrands. Because of the overall factor of u in the integrands in (41 1 it is the terms with 



odd powers of /3 that are now accompanied by an even power of u and so contribute to <Tixy{Ti), whereas the terms even in /3 
vanish after the integration with respect to u. When /3 = 0, <Jixy{Ti) vanishes and aixx{Ti) reduces to the result of Antezza et 
al. IS) for non-moving plates, as one can see by inspection. 



C. Contribution to the Poynting vector 



In the next section we will find that the contribution to the Poynting vector between the plates enters into our method of 
calculating the Casimir-Lifshitz force. The Poynting vector between the plates 



S = — (E X B) 
Mo 



(42) 



can also be obtained from the correlation function (33 1. Similar to (16i and (17i, we obtain an electric-magnetic correlation 
function given by 



A^o Jo 



duj uj ebccN'c (A A) ad (r , r' , t^) ■ 



(43) 



The object of interest is the contribution of Plate 1 to the x-component of the Poynting vector, which we denote by Six{Ti). 
From (42 1 and (43 i we see that Six{Ti) is determined by (33 i as 



Mo Jo 



dui UJ exab^bcd V^(AA)ad(r,r',t^) 



(44) 
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Inserting ( 33 i in (|44|, and again making use of the identities ( 35 1-( 37 1 and the definitions ( 38 i-( 39 1, we obtain 



h 



oo />oo pes 



du dv du) coth 



IGtT^ ./-oo J-oo Jq 

lm{rE2)\aBB\^{cs^ -upuj)^ +lm{rB2)\aBE\^v^P^c^w^ 



2 , ..2 ^2^2^21 



X <^ 4Im(r£;i)- 



2 , (b^s) 



oc poo poo 



CS^ — uPujYaEEdBB + aEBdBEv'^P'^C^w'^ 



„ , du I dv I duj coth I ) uj lics^ - u/Sto)'^ + v'^fi^c 

167r-* ,/-oo J-oo Jcs \2kBTi I 



2 21 
W 



X { [l-\rEi\ 



1 - |rB2|') |aBBp(cs2 - upio)^ + (l - |rB2|') |aBBp?^'/3'c' 
|(cs2 - uPioYaEEaBB + aEBaBEv'^P'^c^w'^f' 



,2^2 



(b ^ b) 



(45) 



Note that the integrands in (45 i differ from those in (41 1 only by a factor of — ^, a fact that will have great significance for the 
Casimir-Lifshitz force. 



IV. CONTRIBUTION OF PLATE 2 



To obtain the contribution of Plate 2 to the Casimir-Lifshitz force we use the results of the previous section to write down its 
contribution to the stress tensor and the Poynting vector in a frame co-moving with Plate 2. In this frame the radiating sources 
in Plate 2 are at rest and Plate 1 is moving in the negative y-direction (speed —(3 in the y-direction). Once we have obtained the 
contributions of Plate 2 in the co-moving frame we will consider how they determine the contributions ct2xx(T2) and a2xy{T2) 
to the stress tensor in the frame of Fig.[T] 



A. Contribution in co-moving frame 



Consider the changes needed to turn the results (40l, (41 1 and (45 i for <7ixx{Ti), (Jixy{Ti) and S'i^(ri) into <J2^^{T2), 
(j'2xy{T2) and S2x{T2), where the prime refers to a frame co-moving with Plate 2. Clearly, the interchange 1 2 is required, 
as regards the reflection coefficients and temperature. Note that the reflection coefficient s ([27| - ( [28] l are evaluated in frames 
co-moving with each plate so we use the same reflection coefficients as in the frame of Fig.fl] Thus the only change required in 
the reflection coefficients and temperature in (40i, (41 1 and (45 1 is the interchange 1^2. The radiating plate (Plate 2) is now at 
X = a with the non-radiating plate at a; = 0, so compared to the last section this entails the following change in the ^-coordinate: 

X —X + a. (46) 



Neither the stress tensor contributions (40i and (41 1, nor the Poynting vector contribution (45 i, are functions of x, but the 
transformation (|46]) affects the tensor and vector components ifTSl : 




'^xa^xb'^ab — O'xxi " xy 

^x ^xa^a 



(^xy — > '^xa'^yb(yiab — —^x 



(47) 

(48) 
(49) 



Thus, in addition to the other changes, a sign change is required in ( |4T] i and (45 i to obtain the contribution of Plate 2, whereas 
this is not the case for (40 1. The moving plate (Plate 1) now has speed —(3 and the integration variables (40i and (41 1 are those 
we would have in the frame co-moving with Plate 2; hence we must put /3 — /3, uj uj' and u u' everywhere except in the 
reflection coefficients. Noting that w is invariant under uj ^ uj', u ^ u' and defining the transformed value of s by 



(50) 
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we can now write down cr2^^{T2), cr2xy{T2) and S2^{T2): 

r du' r dv r dJ coth ( e~^^\^\\w\ UcP + u'aJf + v^lS^^w^] 

J_oo J-oo Jo \2kBT2J 



h 



. ,Re(r^;i)|aBs|^(cs'2+«'/Jw')2 + Re(r-Bi)|a^;B|V/3Vw;2 1 

X < 4Im(r£;2) — ^ ^ 2 \- [E ^ B) > 

[ \{cs''^ + u'l3ujyaEEaBB + aEBaBEv'^l3'^c^w'^\ J 

du' / dv / dio' coth -f^ w [(cs'2 + u'Pcv'f + v^P^c'w^] 

00 J-oo Jcs' \2KbT2J 

I / .\ (1 + l^^il') l«SBp(cs'2 + u'f3u,r + (1 + \rBif) \aEB?v^^^(?w^ 
X Ml - \rE2\^) ^ ^ ^ 



167r3 



|(CS'2 + u'Puj'YaEEaBB + ttEBaBEV'^P'^C^w'^f 

+ {e^b)}, (51) 



167r3 



du' r dv r d^' coth (^^^\ e-^'^l-'lu' [(cs'2 + u'puj'f + v^(3^cV] 

J-oo J-oo Jo \2kBT2J 



[ {{cs'"^ + u'/3u;')'^aEEaBB + aEBaBEv'^P^c'^w^f J 

^ L I -00 £ (^) [^^"^ ^ "'^^'^^ ^ 

(1 - \rEif) lasBlHcs'^ + u'P^'Y + (1 - ksil") |ai5Bpt''/3'c2w;2 



(i-ki.2|') 



|(cs'2 + u'liio'YaEEaBB + aEBiBEv'^P^c^w^l^ 
+ {e^b)}, (52) 



h 



167r3 



du' r dv r duj' coth ( -^^) e-2«l-la;' [(cs" + u'(3oj'f + v^p^cW 

J_oo J-oo Jo \2kBT2J 



X [lIm(rj,,) ^"'^'^^^^'°^^''^^''' + "'^'^'^'+^"'^'"^^^'°^^'''''^'^'"'' +{e^b) 



167r3 



(l-|r-i^2|') 



r dv r dLo' coth ( ^' Kcs" + u'pio'f + v^p^cW] 

-00 Jcs' \2kBT2J 

(1 - IrBll') |aBBp(cs'2 + u'fiuj'f + (l - |rBi|') |ai;B|2z;2/32c2«,2 



|(cs'2 + u'(3uj'YaEEaBB + aEBttBEV'^P'^c'^w'^l 

{e^b)\, (53) 
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The quantities uj and u occur in the reflection coefficients tei and tbi in (51 
both explicitly and also throug h the definitions ([38|)-([38]l. In their appearance in (|5T 



(53 1; we recall that r^i and tbi are present 



of the integration variables to' and u' in accordance with the inverse of the relations ( 29 1: 



([53] ), UJ and u are to be taken as functions 



u) ^ ^{uj' + I3cu'), u^^{u' + I3uj' /c). 



(54) 



B. Contribution in frame of Fig.[T] 

We must now consider how to use the results of the previous subsection to find 0-2x2; (72) and cr2xy{T2), the contribution 
of Plate 2 to the stress tensor components (j2]l and (j3]l in the frame of Fig. [T] Relating tensor components in different inertial 
frames is straightforwardly done by a Lorentz transformation, but here we are dealing with one of two contributions to a tensor 
(the electromagnetic energy-momentum tensor) and when a tensor is decomposed into a sum of two terms the individual terms 
may not be tensors. If one ignores quantum-vacuum effects and considers the purely thermal electromagnetic energy then it is 
clear that the contribution of each plate to the energy-momentum tensor separately constitutes a tensor, since each contribution 
gives the full energy-momentum tensor when the other plate is at zero temperature. In contrast, the quantum- vacuum part of the 
energy-momentum tensor cannot be separated into contributions from the two plates, each of which could, in some arrangement, 
constitute the entire energy-momentum tensor and be separately measured. The contribution of Plate 2 to the purely thermal 
energy-momentum tensor can thus be Lorentz transformed from the co-moving frame to the frame of Fig. [TJ let us first do 
this and then consider the quantum-vacuum part of the energy-momentum tensor Separation of quantum-vacuum from thermal 
effects in the foregoing is achieved by the identity 



coth 



2kBT 



sgn(w) + 2sgn(u;) 



exp 



h\Lu\ 



(55) 



where the first term gives the quantum-vacuum part and the second term, containing the Planck spectrum, gives the thermal- 
radiation part |3 1. We denote the thermal part of a quantity by the subscript th and proceed to Lorentz transform the thermal part 

<^2..{T2) and a[^^^y{T2) of ^ and ^ 

jsTtne 



to obtain at\,2xx{T2) and ath2xy{T2). 



In the gap between the platesTtne energy-momentum tensor T^" is given by ifTSll 



1 n2\ 



S/c 



S/c 
cr 



(56) 



where S is the Poynting vector (42 1 and cr is the Maxwell stress tensor ([T]). (The conventional definition ([T]l takes the Maxwell 
stress tensor as minus the spatial part of the energy-momentum tensor [ 13 1.) We have been writing indices in the lower position 



on S and cr, which is consistent with (56 1 as we choose a metric signature of +2: 



g^u = diag(-l, 1,1,1). 



(57) 



Equation (56i is written without primes and refers to the frame of Fig. [T] In the frame co-moving with Plate 2 the energy- 
momentum tensor has the value T'^'', related to T'"' by ifTSll 



7 7/3 

10 

7/3 7 

1 



(58) 



From (56 1 and (58 1, we see that cr^j^jixC^a) and cr[]-,2xyi'^^) ^re related to c7th2xx{T2) and ath2xy{T2) by 

(Jth2xx{T2) = Cr'th2xx{T2), Cr2xy{T2) = 7 Wth2xy{T2) - /^^thx (^2 ) /c] . 



(59) 



Thus no change is required in the thermal part of (51 1 to obtain crth2a;2;(72)- In contrast, (Jth2iy(72) is given by the linear 
combination (|59| of the thermal parts of (|52]) and (|53|l. This linear combination has a very simple effect, however, because the 
integrands in pll and (53 i differ only by a factor of — we see that (59 1 gives a factor in the integrands of 



7(u' + I3uj'/c) = u. 



(60) 



Thus the only change required in the thermal part of (52i to obtain crth22;a(72) is to replace the overall factor of u' in the 
integrands by u. 
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When adding the contributions from each plate to find the stress tensor components Q and ^ we should use the same 
dummy integration variables in each contribution. We therefore also transform the integration variables ui' , u' in the contribution 
of Plate 2 to w, M using the fact that the Lorentz transformation (29 1 has unit Jacobian ifTSl . The factor cs'^ + u' /3uj' that appears 



in (51 1-(53 1 can be written in terms of unprimed quantities through the identity 

„'2 



cs ■ 



u(3u!. 



(61) 



We must also consider the limits on the frequency integral when we transform in (51 1 and (52 1 to uj. The Lorentz transfor- 
mation in question does not change the value of w and the limits cs' in the integrations over uj' , which represents the separation 
between real and imaginary w, become a limit cs in the integrations over uo. From (29 1, the limit in the integrations over oj' 

corresponds to a limit /3cu in the integrations over u. The integrations duj' over positive values of uj' in the evanescent- wave 
parts of (jsTj) and ( 52 1 thus turn into integrations J^^^^ duo which range over positive and negative values of w, since u takes values 
from —CO to oo. Overall, performing the change in integration variables, (59 1 and (61 1 give (Jt\i2xx{T2) and at\i2xy{T2) as 



{T2 



167r3 



du 



dv dio2 

— oo J f3cu 

2 I „,2o2 2, 



exp 



41m(r£2) 




Re(ri=;i)|aBB|2(cs2-M/3w)2. 



hiu' 

fcs72 



Re(rBi)|a£B|'w'/3'c2«;2 



|(cs2 - ufiujYaEECiBB + QEBaBEv'^P'^c^w'^l 



+ {e ^ b) 



du 



1 - \rE2\ 



dv 



dLu2 



exp 



fcs72 



r/ 2 d \2 , 2 a2 2 2 

w [{cs — UfjuJ) + V p C W 



1 + \rEif) \aBB?{cs^ - upujf + (l + \rBi\^) \aEB?v^ fi^c'^w^ 

|(cs2 - u/3u!)'^aEEaBB + aEBClBEv'^fS'^c'^W^f 



+ {e^b) }, 



(62) 



0'th2xy 



167r3 



oo poo pes 

du dv du!2 

00 <y —00 J j3cu 

.2 a2 2„,,2" 



exp 



kBT2 



X [(cs^ - ul3ujf + v^/3^c^w^ 



4lm{rE2) 



2„,2 «2„2„,,2 



lm(rj;i)|QBB|2(cs2 - uf3uj)'^ + lm(rBi)|Qi;Bp«^/3 
|(cs2 - ufiLoYaEEa-BB + aEBO'BEV^P^c'^w^l'^ 



{e ^ b) 



h 

167r3 



du i dv I du;2 

-00 J cs 



exp 



kBT2 



[l-\rE2[ 



u [(cs2 - u(3ujf + v'^P^c^w^ 

(1 - Vei?) \aBB?{cs^ - upLof + (1 - Irsil') \aEB?v^ P^c^w^ 
|(cs2 - uPwYaEEdBB + aEBaBEv"^ f3'^c'^w'^f 



(e ^ b) 



(63) 



We can separate the integration over negative values of to in the evanescent- wave parts of ( 62 1 and ( 63 1 through 



du I duj = 

1 J f3cu 



^.0 

du I duj - 

Pc\u\ 



du / du) 
Jo 



du / dw, 

Jl3c\u\ 



(64) 



where we assume that /3 > 0; it is easy to verify that a similar treatment of the case /3 < gives the same final result (67 1 
and ([68]l. Through a variable transformation uj —oj, u —u in the first term on the right-hand side of (64i we can write 
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the results ( 62 1 and ( 63 i in a form where all integrations are over positive frequency. To do this we need to consider the effect 
of w — > ~uj,u —u on the integrands in the evanescent-wave parts. Note first from ( [29] l that this transformation implies 
uj' —uj'. We thus need to consider changing the sign of the frequency at which the reflection coefficients of both plates are 
evaluated. A change in sign of the frequency at which they are evaluated is equivalent to taking the complex conjugate of the 
reflection coefficients, whereby the real parts are unaltered but the imaginary parts change sign. From this one can see that a 
simultaneous change in sign of lo, u and lo' in the integrands of the evanescent-wave parts of ( |62| and (63 i produces an overall 
change in sign of the integrands together with a factor of —1 in the exponent of the Planck distribution factor. After the change 
— > — w, u —u, the first and last terms in ( 64 1 can thus be combined to give an integral du J^'^ duj whereby ( 62 1 and ( 63 1 
can be written 



167r3 



du I dv I dijj2 sgn(a;') 

-oo Jo 

..2 q2 2„,.2- 



exp 



h\uj'\ 

kBT2 



- 1 



X i 41m(r£;2) — ^ '-^ ^ 2 h (b ^ -b) 

y \{cs'^ - uPioYaEEaBB + aEBaBEv'^P'^c^w'^] 

-1 



167r3 



du i dv 



— OQ J CS 



dLu2 



exp 



fcsT2 



1 



1 - \rE2\ 



(l + Vei?) WbbWcs^ - ufiujf + (l + \vBif) \aEB?v^P^c^w^ 



(cs^ - u(3wYaEEaBB + asBaBEv'^P'^c'^w'^^ 



+ {e ^ b) 



(65) 



nOO nOO nCS 

Ot\i2xyiT2) = -7—^ I du I dv I (iw2sgn(tj') 

iDT^ J-oo J-oo Jo 



exp 



\kBT2 




lm(rgi)|QBgp(cs^ - uI3ujY + lra{rBi)\aEB?v'^P 
|(cs2 - uPuYaEsaBB + aEBO.BEv'^P^c'^W^ 

-\ -1 



2„,2 fl2„2„,,2 



(i-kiJ2r 



(e ^ b) 



du I dv dui2 

oo J cs 



exp 



kBT2 



1 



u [{cs^ ~ uf3ujf + v'^P'^c^w^^ 



(l - \rEi\') \aBB?{cs^ ~ uPuf + (l - |rBi|') \aEB?v^ P^c'w^ 
|(cs2 - uf3uj)^aEEaBB + aEBO-BE'"'^ P'^c?w'^\^ 



+ (b <-> b) 



(66) 



The thermal parts of the components (|2]) and ([3]l of the stress tensor between the plates are determined by ( [65] l, (|66]l and the 
thermal parts of (40l and (41 1. This is the effect of real radiation reflecting back and forth in the gap between the plates, though 
because of its motion the radiation in the gap arising from Plate 2 does not have a thermal distribution in the frame of Fig. [T] 
The plates have a similar effect on the quantum-vacuum radiation in the gap, but such radiation has a featureless "spectrum" in 
that all frequencies are present in their zero-point guise. If the real radiation is removed then, in the frame of Fig.[T] the effect of 
zero-point radiation is obtained by replacing the spectral distribution of the real radiation from each plate (with an extra factor 
of 2) by sgn(tj), as is seen from (55 i. From (65 i and (66 1 we thus obtain the full contributions in the frame of Fig.[r| a2xx{T2) 
and (T2xy{T2), as 
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0'2xx{T2) 



167r3 



/oo poo pes 

du dv duj <1 -\- 2sgn(cL;') 
-oo J — oo J 



X e-2"l-l \w\ [(cs2 - + v^p^c^w^] 

X I Alm{rE2) — ^ '-^ ^ 2 h (b <^ b) 

[ \{cs'^ ~ u(3LoYaEEaBB + aEBdBEv'^P'^c^w'^] 

n -1 ~1 



oo poo 



oo J —oo J cs 



du I dv I dw {1 + 2 



167r3 

(l + \tei\^) \aBB?{cs^ - u(3u:f + (l + \rBi\^) \aEB?v^ P^c^w^ 



{l-\rE2i 



{cs^ — u(3ujYaEE0-BB + aEBO-BEy^P^C^W^^ 



(e^b)}, 



(67) 



h f°° f°° I 

(^2xy{T2) ^ J du J dv J < 1 + 2sgn(cj') 



exp 



h\Lu'\ 

kBT2 



- 1 



lm{rEi)\aBB?{cs^ -u(3ujf +lm{rBi)\aEB?v'^0^c^'-'^ 
|(cs2 - u(3uj)^aEEaBB + aEBO.BEv'^P^c^w^ 



X |4Im(rB2) 



{e ^ s) 



du / / du) ■11 + 2 

oo J cs I 



exp 



- 1 



167r3 

X u [(cs^ - u/3a;)2 + v^^'^c^vu'^] 

(l - Irsil') |aBBp(cs2 - uPlo)^ + (l - {rsif) \aEB\'v^ P\^w^ 



1 - \rE2\ 



|(cs2 - ufiuYaEEaBB + aEBaBEv"^ P'^c^w'^ 



{e^ b) ). 



(68) 



The components of the stress tensor between the plates that give rise to the Casimir-Lifshitz force are now found from Q and 



V. THE CASIMIR-LIFSHITZ FORCE 

The electromagnetic force on a section of material is given by differences in the electromagnetic stress tensor on the boundaries 
of that section. To obtain the force we therefore require the stress tensor in the frame in which the medium is at rest. The results 
of the previous section were aimed at finding the stress tensor in the rest frame of Plate 1, and without loss of generality we 
consider the force on this plate. 

The calculations in this paper were performed for infinite half-spaces, as in Lifshitz's pioneering work |[l] O. Lifshitz's 
formalism gives the stress tensor due solely to the presence of the cavity, which vanishes when the plates are infinitely far apart 
(a oo); his results are valid for plates that are thick enough for the cavity boundary to have a negligible effect on the local 
electromagnetic fields at the external surfaces. A consequence of the different method of calculation used here and in [6 | is 
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that contributions to the stress tensor that are independent of the plate separation a are included, even in the case considered 
by Lifshitz (non-moving plates at the same temperature), and these contributions are of three kinds. Firstly, we have not yet 
performed the regularization necessary in all calculations of the Casimir-Lifshitz force and so the a-independent diverging stress 
present in the absence of the plates must be removed; for the stress tensor components relevant for the force this divergence 
occurs only in the propagating-wave contribution to axx and it will be isolated in the next subsection. Secondly, the formalism 
used here takes account of the a-independent contribution to the stress tensor that each infinite half-space would make in the 
absence of the other. On its own, each infinite half-space in thermal equilibrium would experience a radiation pressure on its free 
surface, and this pressure, equal to ■K'^{kBTi)^ /{A^1t?c') for Plate 1, is present in the contribution to a^-x as computed here |6|. 
For non-moving plates (/3 — 0) at the same temperature (Ti = T2) the (regularized) result for axx obtained here and in |6| differs 
from that of Lifshitz by these radiation-pressure terms (see ^ for a detailed discussion). Thirdly, in the situation considered 
here and in |6 1, the heat transfer between the plates causes a momentum transfer that is independent of the plate separation, and 
this momentum transfer also shows up in the stress tensor |6|. 

If the plates are taken to have a large but finite thickness the total force requires consideration of the stress tensor on the 
external surfaces, which depends on the details of any thermal radiation in the external region |6|; for example, radiation in the 
external region that is in thermal equilibrium with the plate will cancel the radiation pressure in the cavity due to that plate. 
For the moving plates considered here, one may also need to consider lateral motion of the plate relative to external thermal 
radiation. Our concern here is with the force arising from the stress tensor in the cavity, as this gives the Casimir-Lifshitz force. 

Before we proceed to discuss the force on Plate 1 due to the stress tensor in the cavity, let us consider briefly the electromag- 
netic stress tensor inside the plate. There are always non-vanishing components of the stress tensor inside the plate, even at zero 
temperature, namely (jyy and a^z, and the zero-temperature electromagnetic energy density is also non-zero in the material. The 
non-vanishing zero-temperature stresses and energy density decrease to zero inside the plate as one moves infinitely far from the 
boundary, since the the local electromagnetic field infinitely far away is unaffected by the boundary. The components ayy and 
(Jzz do not cause any forces, but the energy density in the plate contributes to the total Casimir-Lifshitz energy, the derivative of 
which with respect to the plate separation a is minus the perpendicular force. In |5| it was shown that at zero temperature the 
components of the stress tensor relevant for the forces (the xx- and xy-components) vanish inside the plate. The presence of 
thermal radiation, however, produces a non-zero xy-component of the stress tensor in the plate as well as in the gap. To see this, 
note from ([T), Q, ^ and 



D(r,^) = eoeME(r,c^), B{v,lo) = ^lo^l{u;)f^.{v,Lo), (69) 



that Uxy can be written 



axy = {EyDx + BxHy), (70) 
including inside the media of the plate. The operator version of the basic continuity properties of electromagnetic fields ifTSll 



shows that (70i is continuous across the plate boundary. The presence of a non-vanishing, constant axy between the plates 
therefore implies a non-vanishing axy inside the plate, and in both regions the non-vanishing axy is due solely to the thermal 
radiation. If Plate 1 has a finite thickness and is not in motion relative to any thermal radiation in the external region, then the 
component axy inside the plate is due solely to incoming thermal radiation from the gap between the plates. The radiation from 
the gap is absorbed as it moves into Plate 1 and we assume the plate is thick enough for axy to be effectively zero on the external 
surface. The lateral force per urut area on a slice of Plate 1 between x = xi and x — X2, say, is axy{xi) — axy{x2) and since 
axy is effectively zero on the external surface there is a lateral force per unit area in the y-direction given by the (constant) value 
of axy in the gap. The lateral force must of course oppose the relative motion and so act in opposite directions on the two plates. 
To obtain the force on Plate 2 one should consider the stress tensor in its rest frame, but the fact that the lateral force acts in 
opposite directions on the two plates can be seen from axy being constant in the gap since it represents the flow of j/-momentum 
in the (positive) x-direction 1 15 1, i.e. the direction out of Plate 1 but into Plate 2. 

Consider now the expressions for the perpendicular and lateral forces on Plate 1 due to the electromagnetic stresses in the 
cavity. Substitution of ( |40)i and (67 1 into ^ gives the perpendicular pressure axx{Ti,T2); no great simplification results from 



the addition of ( [40| and (51 1 so we do not explicitly write the expression for axx{Ti,T2). The lateral force per unit area is given 



by ([3]), ( [4T] ) and (68 1; because of the symmetry of the integrands in (41 1 and ( [68] ) under the interchange 1^2 everywhere except 
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in the first factors, we can write this as 

<^xy{Ti,T2) 



h 



OC POO PCS I 

du dv diu < sgn(u;') 

-oo -/ — oo •/ I 



xe-2-l-lu[(cs2-u/3a;)2 + i;2/32^: 



exp 



2 21 

■ W 



h\uj'\ 
kBT2 



- 1 



2„,2 o2„2„,,2 



] Ai I Jm(r£;2)|asBp(cs2 -'u/3w)2 + Im(rB2)|asispw^/3^ 

X < \\m[rE\) Y 

I |(cs2 - uiibjYaEEaBB + aEBaBEv'^(3'^c^w'^\ 



1 -1 



{e ^ b) 



+ 



OC />C30 /"OO 

-oo J —QQ J CS 



exp 



^572 



exp 



□2 2 21 



87r3 

(l - \rE2\^) WbbIHcs^ - uPlo)^ + (l - |rB2|') WbeI'v'P'- 



X < 1 - r^si 



(CS2 - uPwYttEEaBB + ttEBttBEV^ 



+ {e ^ b) 



(71) 



The perpendicular force o'2:2;(T'i, 72) consists of a quantum- vacuum contribution plus a contribution from the thermal radiation 
of the plates, as is seen from (40 1, (55 1 and (67 i. The quantum-vacuum part of <Jxx{Ti,T2), which gives the exact perpendicular 



pressure at zero temperature, has a fairly simple expression |5| that will be derived again in the next subsection. Note that the 
lateral force ( 7 1 1 is purely an effect of the thermal radiation and vanishes at zero temperature, as shown in 1 5 1 . Thermal radiation 
provides the means of restoring thermodynamic equilibrium, which requires the damping of the relative motion through a 
lateral force, even when the plates have the same (non-zero) temperature. We see from (71 1 and that the non-vanishing of the 



lateral force when the plates have the same (non-zero) temperature is due to the Planck distribution associated with each plate 
appearing as a function of the co-moving frequency. The vanishing of the lateral force at zero temperature can thus be viewed 
as a consequence of the Lorentz invariance of quantum zero-point radiation — it has the same "spectrum" in every inertial frame. 



A. The quantum-vacuum contribution 



As shown above, the quantum-vacuum contribution to the Casimir-Lifshitz force is directed perpendicular to the plates and is 
given by (Txx{0, 0); from (jlj), (jioj), (67 1 and (55 1 this can be written 



h 

47r3 



du 



dv / dt^e-2'^l'"l|w| 



|(cs2 - uPLoYaEEdBB + aEBO-BEv'^ P'^C^w'^f' 
2„,2 fl2 2„,,2 



X {lm{rEirE2)\aBB\'^{cs^ - ul3ujY+lm{rEirB2)\aBE\^v'^l3'^c^w'^ + [e ^ b)) 



du 



dv 



dujw- 



|(cs2 - uPuj)^aEEaBB + aEBaBEV^f3^c'^W^\ 
X { (l - \rEirE2\^) \aBB\^ics^ - u/Stof + (l - Ir^irssl') WbeI^v^ P^cV + {e 



(72) 



Recall that the two (triple) integrals in ( 72 1 are distinguished by w being imaginary in the first but real in the second. This fact, 
and the overall factor of w in each of the integrals, allows us to write (72 1 as the real part of one (triple) integral. We do this by 
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means of the identities 



Re (e2™V£iri;2aB£;) + \ Wee? = \[^- VEirE2?) (unreal), 
Re (e^^VBirsaa^js) + ^ lassl' = ^ (l " \rEirB2\') {w real), 

Im {e^''"''rEirE2a*EE) = e"2|-l'^Im (rBir^s) imaginaiy), 
Im (e2™Vi5irs2aBB) = e-^l-l^im (r^irsa) (w^ imaginary), 



(73) 

(74) 

(75) 
(76) 



which show that ( [72| l is 

CT..(0,0) 



47]-3 



Re 



/ (iw / duj 



CO J~oo Jo \{cs'^ — uPuoYaEEaBB + a-EBO-BEv'^P^c^w'^Y 



e~'^''rEirE2a*EE + ^ laEE^ ] laBB^ics^ " ^Pujf 



e'^VBirssa^jB + \ \aEB\' ] \aBE?v^ P^c^w^ 



+ {e b) 



(77) 



By multiplying the numerator in (77 i into the factor in braces, we find terms proportional to (cs — u[3uj) v [3 c w , these can 



be re-written using the following (admittedly non-trivial) identity: 

e'''"''rEirE2aEE + \ Wee?^ Wbb]^ + (^e^^VsirsaaBs + \ \aEB\'^ |as£;|' + (s <-> s) 

■ e^™° {rEl1'E2aBBa*EB'^*BE + '>'Elf'B20.BEO*EEO-*BB) + O-EEO-BBaEB^BE + 0*EEO*BB^EBa-BE 

+ {e ^ b). 



(78) 



Using (78 1 in the terms proportional to (cs — u(3ui) v (3 c w in (77 1, it can be written 



f7x.(0,0) 



Re / du I dv I duj 



47]-3 



-cc- JO 



|(cs2 - ufiujYaEEaBB + aEBaBEV^P^c'^w^\ 
X je^*""" [(cs^ - u(3u:Y{rBirB2aEE + rEirE2aBB) + {rBirE2a,EB + rEirB2aBE)v'^ P'^c^w'^] 



,1/2 ,3\2 I 2a2 22|2 

+ (cs — upcj) aEEO-BB + aEBCLBEV p c w \ 



(79) 



which immediately reduces to 



CT2;x(0,0) 



47r3 



Re 



du 



dv 



duj 



(cs^ — u(3ijjYaEEaBB + aEBaBEv'^P^c'^w^ 



X [(cs^ - u(3ujY{rBirB2aEE + rEirE2aBB) + {rBirE2aEB + ?'£;i?'B2aB£;)i'^/3^c^w^] 



47^3 



Re 



/•OO 

dw / duw. 
Jo 



(80) 



The second integral in ([80| is the diverging zero-point stress in the absence of the plates and must be dropped to obtain the 



Casimir-Lifshitz force |3|. The quantum-vacuum force is thus given by the first integral in (80i, which reproduces the result 



of |5|. As in |5 |, we can write (80i as an an integral over imaginary frequencies, where the integrand is better behaved. The 
permittivities and permeabilities are analytic on the upper-half complex-w plane | |T4) and w acquires a positive imaginary part 



in this region, so the integrand in the force term in ( 80 1 is also analytic on the upper-half complex-o; plane. The uj integration 



contour from to cx) can therefore be stretched into the upper-half complex plane so that it runs from to ioo and then back 
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to oo on the real axis in a quarter circle. Since the integrand has no contribution from infinite frequency (where the materials 
are transparent), the quarter circle in the contour makes no contribution and the integration can be written from to ioo. On the 
positive imaginary frequency axis the permittivities and permeabilities are real [14,1 and w is (positive) imaginary; the integrand 
is therefore real and we can remove the Re in ( 80 1. Introducing the notation 



-1 -1 2a|tu| 

r f> I I 

' El' E2^ 



Abe 



i-EB — ' El' B2 

the force term in ( |80l ) then takes the form obtained in ||5l: 

CT..(0,0) 

'{Abe 



CO 

c 

1, 
1, 



IK, 

Abb 

^BE 



-1 -1 2a|ji>| 
' Bl' B2^ 



^BR — ' Bl' E2^ 



1, 
1, 



he 

47^3 



du 



dv \w\ 



+ Abb){s'^ - iK-upf - {Aeb + Abe)w'^v'^I3'^ 



AEEABBis"^ - inul3Y - AebAbew^v^P^ 



(81) 

(82) 
(83) 



(84) 



VI. CONCLUSIONS 



We have obtained the exact solution for the Casimir-Lifshitz force on parallel plates with arbitrary constant lateral motion, 
where the plates are at different temperatures and each has arbitrary electric permittivity and magnetic permeability. This 
generalizes the zero-temperature solution 131 and the solution for non-moving parallel plates at different temperatures The 
perpendicular force has contributions from the quantum vacuum and from the thermal radiation. There is also a lateral component 
of the Casimir-Lifshitz force unless both plates are at zero temperature. The formula for the purely thermal lateral force indicates 
that the non-existence of a quantum-vacuum friction between parallel plates is to be attributed to the Lorentz invariance of 
quantum-vacuum energy. 
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APPENDIX A: EXPANDED FORM OF THE GREEN TENSOR 



The Green tensor given by ( 23 1-( 32 1 contains the inverse matrix 

(l-e^^-^RiRa)"^ (Al) 
which must be expanded before computing the stress tensor. The method for doing this, which includes several consistency 



checks, is detailed in |5 1. The inverse matrix ( Al i was expanded in |5| using a basis constructed from the vectors nE2 and 11^2; 
here it is more convenient to use a basis constructed from n^i and n^i, as follows: 



1 



2iwa 



R1R2) =1 + M, 



M = cee n^i (8) n^;! + cbb (g) + ceb (g) + cbe n^i, 



(A2) 
(A3) 



The unknown c-coefficients in ( A3 1 are found by taking successive dot products of ( |A1| | with polarization vectors; for example, 
Cbb is given by 



Ceb 



R'iR'2 



ab 



(n+Jb(n£;i)a. 



Defining A, i' and p by 



A = n^i • nt 



cs^ — u/Suj 



E2 



(A4) 



(A5) 



n^2 = "£^1 • "52 = - (A6) 
p = [1 + e*™"r£;irB2?-BirB2 - e2™°(r£;ir£;2A2 + rBirB2>-^ + rE2rBi'^^ + rEirB2i'^)] (A7) 



^Bl 
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we write the results for the c-coefficients: 

CEE = pe2"™ri5i (-e2"™r£;2rBirB2 + rE2>? + rB2v'^) 

„2iwa„ „ „ I „ \2 I „ ,,2\ 



~e--'^rB2rEirE2 + rB2^ + rE2V ) , 
CBE = /Oe2™°rsi(r£;2 - rB2)>^v, ceb = P e'^™''rEi{rE2 - rB2)Ai/. 



From ( [23| ) we see that the Green tensor also contains 

R2 (l-e2™"RiR2)"\ 
Using the above results for ( |A1[ ), this can be written in the form 



R2 (1 - e^^'^'^RiRa) 

= d-EE nsi (g) n_Ei + dBB n^i ^ n^i + dEB n_Bi ® n;^^ + dBE n^^ (g) n_Ei, 



where the d-coefficients are 



dEE = e-^'^-CEE/rEi, 
dBE = e-^'^''cBE/rBi, 



dBB = e-^'^^CBB/rBi, 
dEB = e-^''"-'cEB/rEi- 



(A8) 
(A9) 
(AlO) 



(All) 



(A12) 



(A13) 
(A14) 



We insert ( A2 1 and ( A12 1 in (23 1; using the fact that the vector kj^ in (see (20 1) is orthogonal to n^i and nsmi in T (see 
54I1), we obtain 



G(a;, x", u, v, uj) — 



^^'(T + MT) + g-™(^-2a)-iiui2;']>^rji 



2wi 



(A15) 



The Green tensor in the form ( |A15| l was used to calculate the stress tensor and Poynting vector, as described in the main text. 
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